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ABSTRACT 

We present a new method for incorporating nonthermal pressure from bulk 
motions of gas into an analytic model of the intracluster medium in clusters of 
galaxies, which is based on a polytropic equation of state and hydrostatic equi- 
librium inside gravitational potential wells drawn from cosmological dark matter 
simulations. The pressure is allowed to have thermal and nonthermal compo- 
nents with different radial distributions; the overall level of nonthermal support 
is based on the dynamical state of the halo, such that it is lower in more relaxed 
clusters. This level is normalized by comparison to pressure profiles derived from 
X-ray observations, and to a high resolution hydrodynamical simulation. The 
nonthermal pressure fraction measured at r 50 o is typically in the range 10-20%, 
increasing with cluster mass and with redshift. The resulting model cluster prop- 
erties are in accord with Sunyaev-Zel'dovich (SZ) effect observations of clusters. 
Inclusion of nonthermal pressure reduces the expected angular power spectrum 
of SZ fluctuations in the microwave sky by 24%. 

Subject headings: cosmology: theory - large-scale-structure of universe - galaxies: 
clusters: general - intergalactic medium — X-rays:galaxies:clusters - methods: 
numerical 

1. Introduction 



The millimeter and microwave sky has recently been mapped on small angular scales to 
much highe r precision than was previously possible, most notably by the South Pole Tele- 
scope (SPT lLueker et al.ll2010l : IShirokoff et al.ll2011t iReichardt et al.ll201fl ) and the Atacama 
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Cosmology Telescope (ACT iFowler et al.ll2010l ; iDas et al.ll201ll ). These telescopes, as well 
as the PLANCK satellite, will soon provide further results with increased frequency and sky 
coverage. At these scales (multipoles I > 3000) the primordial cosmic microwave background 
is strongly damped, and other components dominate. After accounting for point sources, 
such as dusty star fo rming galaxies and radio galaxies (IHall et al.l 1201 Or IVieira et al.l |2010| ; 
Marriage et al.ll201ll ). the main contributor to microwave anisotropies is then the Sunyaev- 
Zel'dovich (SZ) effect from galaxy clusters, groups, and the intergalac tic medium. SZ mea- 
surements can give strong constraint s on cosmological parameters (e.g. iKomatsu et al.l 12011 



Das et aUbOlll iDunklev et al.ll2010l ; ISehgal et al 



2011: Reichardt et al. 



201ll ). particularly 



the normalization of the matter power spectrum, conventionally given as og, the amplitude 
of fluctuations on an 8/i _1 Mpc scale. 

However, these constraints rely on the ability to predict the SZ signal for a given set of 
cosmological parameters. This requires knowledge of the state of the intracluster medium 
(ICM) in groups and clusters out to high redshift. Many physical processes shape the 
ICM, including shocks from accretion and merging, star formation, supernovae, and active 
galactic nuclei (AGN), among others. Our understanding of these processes and their effect 
on the ICM over time is still evolving. Measurements of how the relations among X-ray 
and SZ observables (and their relation to mass) evolve with redshift are ongoing; curren t 
observational constraints at higher z leave much room for improvement (IReichert et al.ll201ll ). 
On the theoretical side, lacking analytic solutions for gravitational collapse and many aspects 
of the behavior of baryons, recourse is taken to numerical simulations. Full hydrodynamic 
simulations with all the required physical inputs are expensive to run and i nvolve many 



detai l s of subgrid physics; different numerical s chemes can give va rying results ( lAgertz et al. 



20071 : Mitchell et al.l 120091 : IVazza et al.ll2011bl : ISiiacki et al.l 120111 ). Thus, global predictions 
of the expected SZ signal are currently highly uncertain. 

With this in mind, we have deve loped a method for determining the ICM dis t ribution in 



a giv en gravitational potential well ( lOstriker et al.ll2005l ; iBode et al.l 120071 . l2009t iTrac et al. 



201ll ). This makes it possible to take the output of a large volume cosmological simulation 
with only dark matter (DM) and, by populating the halos with gas, produce predictions for 
the SZ signal over a large area of the sky. To date, this model has assumed that the gas is in 
hydrostatic equilibrium (HE) and that it has a polytropic equation of state. Levels of star 
formation and feedback energy are set so t hat the model mat ches current X-ray observations. 
Based on this, the "Standa rd" model of IBode et al.l (120091 . hereafter BOV09) was used to 
produce an all-sky SZ map (ISehgal et al.ll2010h ; this yields a template fo r the angular power 
spectrum which can be compared with mic rowave obser vations. SPT ( Lueker et al. 201ol; 



Shirokoff et al.l l201lt iReichardt et al.l |201ll) and ACT (IFowler et al.l l2010l : IDunklev et al. 



201ol ) both found that the observed power is lower (by up to a factor of 2) than expected 
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from this template. A possible explanation for this discrepancy is that the IBOV09I model 
overestimates the thermal pressure in the outer regions of clusters. X-ray scalings, to which 
the model is normalized, are generally not measured beyond rsoo (the radius within which 
the cluster density is 500 times the critical value) because this region dominates the signal, 
whereas the SZ angular power spectrum is significantly affected by the regions outside this 
radius. 

In particular, bulk motions of gas and turbulent pressure support will occur, making 
the assumption of HE incorrect. This would be more likely in the outer regions of clusters, 
where the gas has been accreted more recently, and in younger, less dynamically relaxed 
clusters. There is considerable support from simulations for t his picture , with the temper- 



ature dropping below that which would be pre dicted by HE (ILau et al.l 120091 ; iBurns et al. 



20101 ; IVazza et al.l l2011al ; iBattaglia et al.ll2011al ). Though difficult to obtain, there is some 
observational evidence that nonthermal sources of pressure are present in cluster outskirts. 
Drops in the average temperat ure with radius, leading to a flattening in the radial pro 



file of entropy, have been seen (Mahd avi et al. 



Hoshino et al. 2010; Urban et al 



2011 



2008 



George et al. 



Akamatsu et al. 



200 



9: 



deviations from sy mmetry (Eckert et al.l 1201 lb 



2011 



2009 



Bautz et al. 2009 



Morandi et a 



Ettori fc Balestra 



this is not always the case llReiprich et al. 
20 111 ; iHumphrey et al.l 120121) . Gas clumping (iSimionescu et al 



2009; 



2011 



201lh. although 



Simionescu et al. 



Urban et al 



2011) 



Walker et al.ll2012l ). and the fact that back- 



grounds dominate (jEckert et al.ll2011al ) complicate interpretation of these observations. 



Trac et al.l (1201 ll ) included nonthermal pressure in the polytropic model, and found 



the resulting SZ template significantly reduced the tension with current SZ measurements. 
However, this had the limitation of assuming that the fraction of nonthermal pressure was 
the same (20%) at all radii, and the same for all cluster masses and redshifts, when in fact 
physical modeling indicates that the nonthermal fraction should be larger in the outer parts 
of clusters, and should also be larger in recently formed clusters (i.e. increase with increasing 
redshift). 

The assumption of a polytropic equation of state means the gas density is given by 

1 r 

p = pod r " 1 and the pressure by P = Pq6 t - 1 , where the polytropic variable 9 is a function 
of the potential, the central density p , and the ce ntral pressure P n (see Eqn.[TJ). Based on 
hydrodynamical simulations of cluster formation, lOstriker et al.l ( 120051 ) assumed P repre- 
sented thermal pr essure only, and argued that the appropriate polytropic index was T = 1.2. 
Shaw et al.l ( 120101 . hereafter SNBL), based on more recent simulations, instead argued that 
it is the sum of thermal and turbulent pressures which follows the polytropic relation (but 
still with T = 1.2). In this case the fraction of nonthermal pressure varies within a cluster. 



increasing with radius, and also increases with redshift. The simulations of IBattaglia et al. 
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( ]2011at ) give support for a cluster mass dependence as well. 



In this paper we modify the m odel of IBOV09I to include a treatme nt of no nthermal 
pressure similar to that presented in ISNBLl . Our approach will differ from ISNBLl in that it 
will be based on the dynamical state of each individual halo, such that less relaxed halos 
have a higher level of turbulent pressure. The relation between the level of virialization and 
the fraction of turbulent pressure is calibrated by comparison to cosmological hydrodynamic 
simulations (in particular a high resolution run using the ENZO code) and the pressure 
profile derived from X-ray observations. 

Sec. [2] summarizes the polytropic model. Sec. |3] presents the new method for including 
turbulent pressure, which includes a free parameter that is set by normalizing to observed 
and simulated clusters. Note the model parameters are set by comparison to X-ray data 
only. The model, applied to a large set of halos, is then compared to SZ observations in 
Sec.H 



2. Polytropic Model 



The goal behind this work is, given a particular cosmology, to produce a simulated 
map of the microwave sky. The strategy is to first run a pure dark matter simulation of a 
large volume, saving the mass distribution along the past light cone. The group and cluster 
sized halos in the light cone are identified, and each in turn is populated with gas; thus the 
thermal and kinetic SZ effect arising from clusters can be calculated. In this paper we will 
use halo sam ples drawn from a light cone output of a simulated cube 1000 /i _1 Mpc on a side 

(0.044, 



sec 



BOV09I for more details) with cosmological parameters Q m , Q\, h, n 



0.264 , 0.736, 0.71, 0.96, 0.80), consistent with the WMAP 7-year results flKomatsu et al. 
201 lh . 



Adding gas to halos is done with the method of IBOV091 Gas is placed in the DM 
gravitational potential <fi such that it is in hydrostatic equilibrium and has a polytropic 
equation of state. The polytropic variable 9 is 







1 + — ^(0o 

r p n 1 



(1) 



where 0o, Po, and Pq are the central values at radius r = 0. IBOV09I assumed that Pq was due 



to thermal pressure alone. Here we will instead follow ISNBLl . so that it is the sum of thermal 
and turbulent pressures which follows the polytropic relation. The turbulent component will 
have the same compressibility as the thermal component, so the polytropic index remains 
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T = 1.2. Thus the total pressure is the sum of thermal P t h and turbulent P nt h components: 



P = P th + P nth = P 9—i 

fnth — Pnth/ (Pth + Pnth) 
Pth = (1 " fnth)P O 0^ 



(2) 
(3) 
(4) 



The density profile is unchanged from IBOV09I . but the temperature is altered by a factor 
(1 — fnth)- fnth varies with radius, so the effective polytropic index, T e ff = d log Pth/d log p 
is not constant, but rather increases with radius. We will let f n th also vary with cluster mass 
and redshift (see Sec.E]). 



Two constraints are used to fix the central values po and Pq (see IBOV091 ). First, the total 
pressure at the outer radius must match the surface pressure P s (calculated from the DM in 
a buffer region at r vir ). Second, energy must be conserved. Assume the gas initially followed 
the DM density distribution inside the virial radius r vir , with the same specific energy. Some 
of the gas (that with the lowest entropy) is removed and placed into stars. The remaining 
gas has an initial energy E g \ this gas is now rearranged to obey the polytropic equation of 
state. The energy of the gas in this final arrangement is Ef — | J PdV (the integral is over 
the cluster volume, out to the outermost extent of the rearranged gas, r^; note Tf can be 
larger th an r vir ). This energy must equal the initial value, modified by various processes (see 
BOVOol for details): 

E f = E g + AE P + e D \E D \+e F M,c 2 . (5) 

The first modification to the initial gas energy is a change due to expansion (or contraction) 
of the gas, AEp = (47r/3)(r^ ir — r^)P s , which accounts for mechanical work done against the 
surface pressure. The next term is a dynamical energy input, needed becau se during cluster 



form ation energy is transferred from DM to gas by dynamical processes ([McCarthy et al. 
20071 ). The binding energy E D of the DM is found, and 5% (en = 0.05) is added to the 



gas; this value provides a good match to the gas fract ions seen in "adiabatic" (no radiativ e 
cooling or star formation) hydrodynamic simulations (ICrain et al.l 120071 : lYoung et al.ll201ll ). 



The final term in Eqn.[5] is feedback from star and black hole formation, where M* is 
the mass of stars inside the virial radius The mass in stars at z = is set to 



MjM mr = A, (5.0 x lO 13 M /M 5O o) ( 



(6) 



where M 50 o and M vir are the tot a l mas s inside r 50 o and r vir , respectively. From the infrared 
properties of clusters, iLin et al.l (120031 ) found A* = 0.026 and a* = 0.26 when meas uring 
the stellar mass inside rsoo- For ne arby, low-mass g alaxy clusters, iBalogh et al.l ( 1201 ll ) find 
stellar fractions in agreement with ILin et al.l ( 120031 ). Note we use th ese parameters to give 
M* for the entire cluster out to the virial radius (not just inside rsoo); lAndreonl (120101 ) found 



-6 - 



that this choice gives a model that matches the observed stellar fract ion within r2n n - Thi s 
stellar fraction differs (particularly at lower masses) from that found by lGiodini et al.l ( 120091 ). 
who obtained A* = 0.05 and a* = 0.37 (again measuring inside r 500 ). However, we find that 



with the iLin et al.l (120031 ) values the model clusters give a better match to observed electron 
pressure profiles (see Sec.|3]). The reds hift evolution of ste l lar m ass, including gas recycling, 
is handled in the manner described in IBOV09I . ILin et al.l (120121 ) find little evolution in the 
stellar fraction relative to ILin et al.l ( 120031 ) out t o z ~ 0.6. At higher re dshifts (~0.2-l), 
and measuring out to 200 times the mean density, iLeauthaud et al. find total stellar 

fractions lower than what we are using here (see also ILeauthaud et al.l 12011 at ); adjusting for 



intra-cluster light would at least partially reduce the difference. 

Feedback energy, from both supernovae and AGN, is parametrized as epM^c 2 . This 
will be set to match observations (Sec. 13.31) ; we find ep = 8 x 10 -6 . If the mass in black 
holes is ~ 10^ 3 M„ this means roughly 6 x 10 3 of the black hole rest mass is deposited 
in the surrounding gas (ignoring the supernova contribution). A similar value of 5 x 10~ 3 
is needed in simulations to produce the correct r elationship between black hole mass and 
stellar velocity dispersion ( iDi Matteo et al.l 120051 ). If 15% of the baryons are turned into 
stars, then this level of feedback adds 1 keV per particle initially inside the virial r adius ; 
this is actually lower than the value of 2.62 ± 0.85 keV per particle IChaudhuri et al.l (120121 ) 
derived from X-ray data. 

It is possible to include relativistic pressure, e.g. cosmic rays and tangled magnetic fields, 



pressure, i.e. P rei = 8 re i{Pth + Pnth) 
will not do so (in effect setting 5 re i 
fit to observed clusters. 



sec 



Bode et al. 


2009; 


Trac et al. 


2011) 



0), because including this effect did not improve the 



3. Incorporating nonthermal pressure 

3.1. Nonthermal pressure from turbulence and bulk motions 

Based on adaptive mesh refinement hyd rodynam ic simulations of cluster formation 
(which include star formation and feedback), ISNBLl found that the radial profile of the 
nonthermal pressure fraction follows a power law in radius: 



nth 



a 



r 500 



n nth 



(7) 
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with n nt h = 0.80 ±0.25. Examining their simulations out to z 
dependent normalization 



1, SNBL derived a redshift- 



ce(z) 



MIN[a (l + zf, (4- n " th ao 1 - 1) tanh(^) + 1] 



where a = 0.18 ± 0.06 and (3 = 0.5. iBattaglia et al.l (1201 laf ) also used simulations to 
constrain f n th, with clusters out to z = 1.5 drawn from a set of smoothed particle hydro- 
dynamics cosmological boxe s; star formation and AGN feedback were included. They saw 
redshift evolution similar to ISNBLl . but also a weak mass dependence; for z < 1 this can be 
included by modifying a to be 



a(z,M) = aoil + zf 



M> 



200 



3 x 10 14 M Q ) 



0.2 



(9) 



The se results seem to be ro bust with regard to the details of the gas physics, as both ISNBL 



and IBattaglia et al.l fcOllaf ) obtain similar values when not including gas cooling and star 
formation. 



This can easily be incorporated into the method of IBOV09I . One possibility is to simply 
use Eqn.[8] or [9] for a. However, it is seen in simulatio ns that relaxed c l usters tend to have 
less turbulent pressure support than unrelaxed ones (lLau et al.l 120091 ; IVazza et al.l 1201 lal : 

201lh . Thus, rather than giving all 



Paul et al. 2011; Hallman k Jeltema 2011; Nelson et al. 



halos at a given mass and redshift the same value for nonthermal fraction, it would be more 
accurate to base a on the dynamical state of the halo. We will do this as follows in the next 
section. 



3.2. Setting the level of nonthermal pressure 

The mean nonthermal energy fraction in the gas is given by / = J P nt hdV/ J PdV, 
with the integrals over the cluster volume out to the virial radius. In a situation where the 
density and temperature distribution of the gas is known (e.g. from an observed cluster or 
a hydrodynamical simulation), the fraction of energy due to turbulent nonthermal pressure 
can be computed directly as 

I y fj,m p 3 J 

where a 2 is the velocity dispersion, k is the Boltzmann constant, and \xm v is the mean mass 
of the ICM particles. 



(10) 
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When determining the polytropic gas distribution for a DM-only halo, the total pressure 
P is found before determining the turbulent and thermal components separately. For a given 
cluster we can thus compute a measure f' ol proportional to the nonthermal energy fraction 
fpoi when assuming the form of Eqn.O 



- _ Ja(r/r 500 ) n ^PdV _ 

J poi j PdV ~~ po1 ' ' 



We wish to base the value of / = af l on the properties of the DM halo; thus a 
measure of how relaxed or disturbed a halo may be is required. The virial theorem states 
1/2 = 2K + W — Eg, where K and W are the kinetic and potential energies, / is the moment 
of inertia, and Es is a surface pressure term (the surface is located at the virial radius r vir ). 
Thus for a DM halo in HE (J = 0) we expect kinetic energy K vt = (Es — W)/2. Let us 
identify any variation from this as being due to an out of equilibrium state, which will also 
be reflected in turbulence. Then we can compute from the DM a measure of the degree to 
which the system is out of equilibrium: 



\K-K. 



vt 



neq 



K 



2K 



E s - W 



(12) 



We will assume that the level of nonthermal support in the gas is proportional to this 
measure, or / = nf neq , where k is a constant independent of mass and redshift. Thus, 
limiting the thermal pressure to be less than the total at the outermost gas radius 77 (i.e. 

a ( r f / r 50o)™ n " 1 < 1), we have 



a 



MIN [ nf neq /f; olJ (r 500 /r f r^ } . (13) 



It only remains to set the constant of proportionality k; once this is known, P n th and Pth 
can be computed in the polytropic model. We will determine k in two ways: by matching 
the observed electron pressure at r 500 , and by computing the thermal pressure directly from 
the gas in a simulated cluster which includes full hydrodynamics. 



Fro m X-ray da t a. lArnaud et al.l (120101 ) derived a electron pressure profile of low redshift 
clusters; ISun et al.l ( 120111 ) found this same profile holds for local galaxy groups. This profile 
is shown as a dashed line in Fig.[TJ Note that, due to the paucity of observations, outside 
of rsoo this profile is derived from simulations; also, for higher redshifts it assumes self- 
similar scaling. To compare the model profiles with this universal form, we divided a sample 
of low redshift (z < 0.2) DM halos from the lightcone into three mass bins: 3 x 10 13 < 
M 5OO /M < 9 x 10 13 ; 9 x 10 13 < M 5OO /M < 3 x 10 14 ; and 3 x 10 14 < M 5OO /M < 9 x 10 14 . 
The mean pressure profile was calculated for each mass bin, and k was varied to minimize 
the difference between the observed and model profiles at r 50 o- For this purpose Pth was 
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converted to electron pressure and normalized by P500 (lArnaud et al.ll2010l ); the observed 
profile was calculated using the mean mass and redshift of the model halos in each bin. The 
best fit value is k = 0.70; the resulting profiles are shown as solid lines in Fig.[T] (the other 
parameters in the polytropic model are as listed in Sec.|2]). 

Steepening the radial dependence, by increasing n nt h to as much as 1.2, did not pr oduce 



a significantly improved match to the universal profile, so we leave n nt h at the [ SNB 



value 



of 0.8. (It may be that n nt h depends on how relaxed a cluster is; see iNelson et al.l 120111 ). 
Cutting ep in half, to 4 x 10~ 6 , also did not change the best fit k. Note the outermost radius 
of the polytropic gas, 77, varies from cluster to cluster. In most cases r/ does not extend 
much beyond 2.6rsoo for the high mass bin; this limit increases to 3.2rsoo for the middle bin, 
while for the lowest mass bin it can extend to 4rsnn- T he gas outside of rj needs to be treated 
separately when making SZ maps JSehgal et al-fcioh . 



As an alternative method of fixing k (and thus the level of nonthermal pressure), we 
compare th e resulting profile to that of a cluster taken from a hydrodynamic al simulation 



( Cenl 120111) performed with an Eulerian adaptive mesh refinem ent code, Enzo (jBryanlll999 



Brvan & Normanlll999l ; lO'Shea et al.l 12004 



Joung et al.ll2009l ). First a low resolution sim- 



ulation with a periodic box of 120 ft _1 Mpc on a side was run, and a region centered on 
a cluster of mass of ~ 2 x 1O 14 M was identified. We then resimulated this subbox with 
high resolution, but embedded in the outer 120/i _1 Mpc box in order to properly take into 
account the large-scale tidal field and the appropriate boundary conditions at the surface of 
the refined region. The subbox centered on the cluster has a size of 21 x 24 x 20/i~ 3 Mpc 3 . 
The initial conditions in the refined region have a mean interparticle separation of 117/i~ 1 kpc 
comoving and a dark matter particle mass of 1.07 x lO 8 /i -1 M . The refined region is sur- 
rounded by two layers of buffer zones, each of ~ l/i -1 Mpc, with particle masses successively 
larger by a factor of 8 for each layer; the outermost buffer then connects with the outer root 
grid, which has a dark matter particle mass 8 3 times that in the refined region. The mesh 
refinement criterion is set such that the resolution is always better than 460/i -1 pc physical, 
corresponding to a maximum mesh refinement level of 11 at z — 0. 

The cosmological parameters of the simulation are Q m , h, n s , cr 8 ) = (0.046, 0.28, 
0.72, 0.70, 0.96, 0.82). The simulation includes a metagalactic UV ba ckground (IHaardt fc Madau 



19961 ) and a model for shielding of UV radiation b y neutral hydrog en (ICen et al. 



2005 . It also 



includes metallicity-dependent radiative cooling ( ICen et al.lll9951 ). S tar particles are create d 
in cells that satisfy a set of criteria for star formation proposed by ICen fc Ostrikerl ( 119921 ); 
star particles typi cally have masses of ~1O 6 M . Supernova feedback from star formation is 
modeled following ICen et al.l (120051 ). 



To emulate the manner in which the polytropic model is normally applied, the DM 
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particles were taken from the halo at selected redshifts, and their masses were adjusted by 
a factor (1 + Oj/Am/) to account for the baryonic mass. These particles were then used as 
input to the polytropic model code. The assumed stellar fraction was increased to match 
that of the simulation, but feedback was not included (i.e. ep = 0) because the simulation 
did not include AGN. This cluster grows in mass from M 500 = 4.1 x 1O 13 M at z = 1 to 
M 50 o = 1.9 x 1O 14 M at z = 0. The thermal pressure profile of the simulated gas at different 
redshifts is shown as dashed lines in Fig.[2j The polytropic gas model without turbulent 
pressure is shown as dotted lines; this gives a larger pressure in the cluster outskirts than is 
seen in the simulated gas. If we instead fix k = 0.7, as determined from the observed pressure 
profile above, the resulting profiles are shown as solid lines; these are much closer to the level 
of the simulation. Alternatively, from the simulated gas one can compute / directly, and thus 
k. Using the value of k taken from the simulated gas gives the dot-dashed curves in Fig.[2J 
At z = (0, 0.5, 1), the simulated gas gives k = (0.57, 0.92, 0.69), for a mean of 0.73 ±0.18. At 
z = 1.6, as the cluster is just forming and is still quite unrelaxed, either choice of parameters 
gives less turbulent pressure than is seen in the simulation. For comparison with Eqn.[7l 
it is also possible to measure f nt h as a function of radius directly from the gas. A power 
law with n nth = 0.8 is a good fit inside r 500 , but in the outer regions substructures cause a 
more variable profile not well represented by a power law. Also, there are asymmetries; if 
we measure f n th in different octants, there is a variation of roughly 30% around the mean at 
the virial radius. 

Based on these two determinations, in what follows we will fix k = 0.7 and n n th = 0.8. 
The effect of basing a on the DM halo state in this manner can be seen by comparing to an 
a wh ich var ies only with redshift. We did this by recomputing the profiles using a values 



from ISNBLl . i.e. using Eqn.[8]with a = 0.18. The resulting profiles are show n as d otted 



as a 



lines in Fig.[TJ We assume a different stellar fraction at lower masses than did ISNBL ; 
result, our mean profile is closer to the universal profile than they found. Basing a on the 
dynamical state of the cluster gives a larger dispersion in thermal pressure profiles with a 
higher mean, particularly at lower masses. However, even with a constant ao the dispersion 
in P t h is large, roughly 50% near ~ 2rsoo, due to variations in DM profiles for halos with the 
same mass. 



For a further comparison to simulations, iBattaglia et al.l (l2011bf ) recently derived a 
fitting function for the mean P t h from a suite of SPH simulations. For low redshifts, this 
mean profile (calculated using the mean M200 and redshift of the halos in each bin, and 
adjusted to electron pressure) is shown as a dot-dashed line in Fig.Q] I n the ou ter regions of 



the halos, this gives profiles similar to those we found when using the ISNBLl form (Eqn.E]) 
for a. Including mass dependence, as in Eqn.[9l does not produce a significant change in 
the mean profile, except in the lowest mass bin where the mean is 11% higher near 2rsoo- 
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lt is counterintuitive th at the mean profile fou nd when using Eqn.[9] is a poorer match 
to the fitting function of iBattaglia et al.l ( 1201 lbl ) than that found when using Eqn.|HJ The 
explanation likely lies in the treatment other physical processes than the nonthermal pressure. 
In particular, the amount of star formation has a significant impact on the entropy of the 
remaining gas. Thus differences in how the stellar mass fraction scales with cluster mass and 
redshift will affect the p and P t h profiles. 

The variation of a with mass at low z (using all halos in the light cone with z < 0.2) is 



shown in Fig.[3J At higher masses, the mean is close to the 



SNBL 



value of 0.18, but there is 



a weak mass dependence in the same sense as seen by IBattaglia et al.l (j2011al ). However, the 
mean a scales roughly as M^q 5 , a slightly shallower scaling than in Eqn.[9j The increase is 
due primarily to an increase in f neq with mass, which is shown in the lower panel of Fig.HJ 
More massive halos, being rarer, will likely have collapsed more recently and thus be more 
out of dynamical equilibrium. We can quantify this expectation by considering the peak 
amplitude l/a(M), where a(M) is the rms density fluctuation on scales corresponding the 
the spherical tophat collapse of halo of mass M. T he formation rate of h alos at a given 



mass and redshift will depend on this parameter (e.g. iSheth fc Tormenlll999[ ). as halos with 
higher peak amplitudes will be more extreme fluctuations. The solid line in the lower panel 
of Fig.H]is directly proportional to the peak amplitude, nf ne g = 0.092/a(M). This tracks the 
mean behavior of the DM halos well, except at the lowest a (highest M). The other factor 
determining a, the change in f' pol with mass, is less important. This change arises because of 
the changing concentration of the underlying DM halos. f ol increases with concentration, 
i.e. with decreasing mass; but as concentration is quite weakly dependent upon M, the 
change is slight (less than 20%). 

The variation of a with redshift is shown in the top panel of Fig.Hl The halos shown are 
restricted to a small mass range so that the evolution shown is unaffected by the change in 
the mass function with redshift. At higher z, clusters are typically dynamically younger and 
less relaxed, hence a is larg er. The increase in a with redshift is slightly more pronounced 
than that found by ISNBLl and IBattaglia et al.l (1201 laj ) ; the mean scales as (1 + z) a76 out 
to z — 3. Again the main driver of this evolution is the variation of f neq , shown in the 
lower panel. Assuming f neq is proportional to peak height reproduces the trend seen out to 
z pd 1.5, but again as a becomes small (i.e. higher z) this leads to an overestimate. Also, 
once again the variation of f' pol is smaller, at about 20% over the whole redshift range shown. 

The evolution with redshift of the mean pressure profile is shown in Fig.[5j A lower 
mass range was chosen so that the mean mass 5.6 x 10 13 M^ does not change with z. The 
model profiles are compared to the lArnaud et al.l (120101 ) profile, whi ch assumes a sel f -similar 
evolution. Our model deviates from this assumption; the fit of IBattaglia et al.l (j2011bl ) 
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deviates in the same manner, but not as strongly. The main processes which can cause a 
break from self-similarity are star formation and feedback; without matching the evolution 
of these processes it is difficult to compare methods. However, we have already seen that 
our method does give a stronger scaling of a with redshift (cf. Fig.HJ). 



3.3. Setting the level of feedback 



The level of feedback, parametrized as ej?M»c 2 , needs to be determined. To constrain ep, 
we fit to X-ray observations in the kT — E(z)M g pla ne, where E 2 ( z ) = Q m (l + z) 3 + • The 



two samp les we compare to are the cluster sample of iMantz et al.l (120101 ) and the lower mass 
sample of ISun et al.l (120091 ) . M g is found inside r 500 , and T is the spectroscopic temperature 
inside the radial range [0.15 — l]r5oo- The simulated sample includes all halos in the light 
cone with M500 > 5.5 x 10 14 M Q and z < 0.45, 79 in all, plus another 79 chosen from z < 0.1 
to span the range 2 x 10 13 < M 500 < 1.5 x 10 14 . 

To find the best fit we minimize y 2 , defined as follows. For simulated halo i and observed 
cluster j (with observational uncertainties <jm and or), 

2 / » . T » r \ 2 

! I -r I ! , (14a) 



a Tj 




(14b) 
(14c) 



The best fit feedback level when including our model for a is ep = 8 x 10~ 6 . If we instead 
use Eqn.[8l the same value is recovered. The X-ray properties of the model are compared to 
the observed samples in Fig.[6j 



This amount of feedback is larger than that in the Standard model of IBOV09I . which 
has no nonthermal pressure. The Standard model value was also chosen to match observed 
clusters in the kT — E(z)M g plane. Starting with this model, the addition of nonthermal 
pressure reduces T without making any change to M g , so at a given T the value of M g is 
now too high. Increasing the amount of feedback both reduces M g and increases T, bringing 
the model back into agreem ent with the observed relation. This amount of feedback is also 
larger than used by ISNBLl . but we are also using a lower amount of star formation. As 
discussed in IBOV09I . the star formation and feedback parameters are somewhat degenerate, 
such that increasing the amount of star formation in lower mass clusters would require a 
lower feedback parameter. Thus a greater amount of feedback is to be expected. 
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As mentioned in Sec.[2l it is possible to also include relativistic pressure, parametrized 
as S re i times the non-relativistic pressure. However, repeating the above procedure while 
varying both ep and 5 re i yields a best fit with 5 re i = 0, so we do not include it. The main 
effect from including 5 re i is to increase the total nonthermal pressure in the core. However, 
there are other, compensating effects we are not including, such as cooling and the presence 
of a cD galaxy. 



4. Implications for the SZ signal 



Having set the parameters of the model by comparison to X-ray observables and sim- 
ulations, it is now possible to make predictions co ncerning the SZ signal. We begin by 
comparing to early results from the Planck satellite. [Planck Collaboration &: et al.l (J2011bJ) 
combine measurements of D\Y = J PthdV inside r^o with X-ray data taken by XMM- 
Newton for 62 z < 0.5 clusters. Y is expected to scale closely with the X-ray quantity 
Y x = CM g T. Here C = g t / \m e c 2 fi e m p ) , where gt is the Thompson cross-section, m e c 2 is 
the electr on rest mass, and ji e m v is the mass p er electron in a fully ionized plasma. The 
data from lPlanck Collaboration &: et al.l (l2011bl ) is shown in the left-hand panel of Fig.[7J the 
gas mass M g is measured inside r 500 , and T is measured in the radial range [0.15-0.75]r 50 o- 
For the f nt h model clusters we include all halos in the light cone with z < 0.5. Fitting 
those with D\Y > 10~ 5 Mpc 2 , the slope of the D\Y — Y x relation is consistent with unity; 
this agrees with iPlanck Collaboration fc et al.l (l2011bT). but is roughly 2-sigma steeper than 



the slope found by iRozo et al.l (120121 ) . If, following iRozo et al.l (120121 ) , we make the pivot 
point 8 x 10~ 5 Mpc 2 , then the model gives \n(D\Y/Y x ) = —0.152 ± 0.002 with rms scatter 
g = 0.031; in other words (D 2 aY/Yx = 0.86, significantly lower than the best fit given by 
Planck Collaboration fc et al.l (l2011bl ). On the other hand, it is higher than the an alysis of a 
subsample of the Planck clusters using different X-ray data by IRozo et al.l (120121 ). who find 
D\Y/Y X = 0.82 ±0.024. These results are insensitive to the redshift range used, but includ- 
ing clusters with D\Y below 10~ 5 Mpc 2 lowers the ratio to (D\Y/Yx = 0.82 and increases 
the scatter to g = 0.037. 

The relation between SZ signal and the X-ray luminosity is shown in the right-hand 
panel of Fig. [71 The luminosity Ly is in the band [0.1-2.4 |keV, measured inside r 50 o- The 
same clusters from IPlanck Collaboration &: et al.l (l2011br) ar e shown, as well as data for 
lower- mass clusters from lPlanck Collaboration fc et al.l (j2011al ). who measured Y for roughly 
1600 X-ray selected clusters and averaged the results in bins of Lx- The model reproduces 
the observed relation well except at the highest Lx, where the model cores are not dense 
enough. Yx depends roughly linearly on gas density, and the core was excluded when com- 
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puting T. In contrast, Lx depends on the square of density, and thus strongly weights 
the cluster center. Our model does not include cooling in the cluster core; such cooling 
would result in higher L x without significantly affecting Y or Y x - The slope of the f nt h 
mod el Lx — Y relation (set mostly by the lower mass clusters) is 1.09, in good agreement 
with lPlanck Collaboration fc et all (j2011al ). For this model, the evolution of Y with redshift 
follows the standard self-similar scaling out to z = 1. 

The above comparisons were all measured within of r^n, whe r eas th e effect of nonther- 
mal pressure is greatest outside this radius. Recently lHand et all (120111 ) used ACT data to 
measure the Y out to a cylindrical radius corresponding to a cluster density of 200 times 
the mean. This includes the entire contribution along the line of sight, as well as a larger 
solid angle than subtended by r 500 , and so is more sensitive to the cluster outskirts. The 
measurement was done by stacking the locations of LRGs, which reside in massive halos, 
to find the projected SZ signal. This data is shown in Fig. [HI Only the radio-quiet LRG 
sample is shown, to avoid contamination. The masses of the halos are determin ed from wea k 
lensing, as the masses determined from bias may be affected by systematics (jMord 1201 ll ). 
The model fits the data within the errors. 



Fig.[8jalso shows the mean relation for the Standard model of lBOV09i . This model has 
no nonthermal pressure and a lower amount of feedback, such that inside rsoo it would closely 
resemble the f n th model. Because of higher thermal pressure outside this radius, it shows a 
larger SZ signal. Including f nth has little effect on the slope of the relation, but reduces Y 
at a given mass by 15% from the Standard model. We also show the nonthermal20 model 
of iTrac et all (120111 ); this has constant 20% nonthermal pressure at all radii (note it has a 
different level of star formation and feedback as well). Thus, while it has a higher thermal 
pressure at large radii, it will have more nonthermal support in the core, and thus the central 
density and thermal pressure will be lower. This leads to a smaller value of Y at a given 
mass. However, the slope of this relation is again little changed. 

The reduced SZ signal predicted when f nt h is included will alter the predicted an gular 
power spectrum in SZ maps. Computing the C\ in the same manner as lTrac et all (120111 ). the 
f n th model yields a thermal SZ template 24% lower than the Standard model (at 280 GHz, 
for cr 8 = 0.8). Two effects which one might expect could alter the shape of the template, the 
generally lower P n th in cluster outskirts and the variation in the level of f n th from cluster 
to cluster, appear not to do so; the peak of the SZ power, at I ~ 3000, remains the same. 
This is in contrast to the nonthermal20 model, which is 45% lower than the Standard model, 
and has a pea k at slightly larger angular scales. The reduction of 24% from the Standard 
model used in Sehgal et al.l (2010) w i ll alleviate somewhat, but not eliminate, the tension 



with observations (jSehgal et al 



20101 : iDunklev et al.ll2010l : iReichardt et al.l 120 111 ). 
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To summarize, in this paper we have pre sented a method of including turbulent pressure 
in a model of the ICM, similar to that of ISNBLl . The fraction of pressure in turbulent 
form is allowed to vary from cluster to cluster, depending on the host halo's dynamical 
state. The normalization is set by comparing the pressure profiles to X-ray observations 
and hydrodynamical simulations. The nonthermal pressure fraction measured at rsoo is 
typically in the range 10-20%, trending higher with cluster mass and with redshift. This 
will be useful in creating improved templates for interpreting new data from millimeter and 
microwave instruments such as the PLANCK satellite. However, a better understanding of 
the thermodynamic state of the ICM at higher redshifts will be needed to fully exploit this 
new data. 
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Fig. 1. — Electron pressure profiles of model clusters. Each panel is labeled by the mass 
range u sed; the mean r edshif t is z ~ 0.1. The dashed line is the universal profile from X-ray 
data by lArnaud et al.l (120101 ) . who give an expression for the normalization factor P^oo- The 
solid lines indicate the mean and one standard deviation for the model clusters when basing 
a on the dynam ics (Eqn.[13]). The dotted lines instead use const ant an = 0.18 in Eqn.p 
following SNBL . Dot-dashed lines are the mean simulated profile of Battaglia et al. ( 2011bf ). 
based on the mean M 2 qo and redshift of each bin. 
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Fig. 2. — Comparison of radial profiles of electron pressure from a hydrodynamical simulation 
(red dashed lines) to the polytropic model applied to the DM potential. Each panel is labeled 
with the redshift. The x marks the virial radius. The dotted line is the model without 
nonthermal pressure. The solid line includes it, with k = 0.7 in Eqn.[13j The dot-dashed 
line instead uses the k parameter measured directly from the simulation at each redshift. 
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Fig. 3. — Top: a = P nth /P at r = r 500 , as a function of mass, for halos with redshift < 0.2. 
Points with error bars sho w the mean and standard deviation in logarithmic mass bins. The 
dotted line is the SNBLl value at z = , an = 0.18 ± 0.06. The dot-dashed line includes 
mass dependence as in iBattaglia et all (j2011al ). Bottom: nf neq of the halos. The line is 
proportional to the peak height l/a(M). 
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Fig. 4. — To;?: a as a function of redshift, for halos in the range 4 x 10 13 < M/M < 8 x 10 13 . 
Points with err or bars show the mean and standard deviation in redshift bins. The upper line 
us es the SNBLl value, an = 0.18±0.06; the lower line adjusts this by the mass-dependent term 
of iBattaglia et al.l (j2011al ). Note that these latter two fits are based on z < 1.5 clusters, and 
the lines are extrapolated beyond this. Bottom: Kf neq of the halos. The line is proportional 
to the peak height l/a(M). 
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Fig. 5. — Electron pressure profiles at different redshift. Solid lines: the mean adopted model 
profile, for selected halos in the mass range 3 x fO 13 < M 5OO /M < 9 x 10 13 and in redshift 
bins 0.4 < z < 0.7, 0.7 < z < 1.0, and 1. < z < 1.3 (panels a re labeled with mean redshift). 
Dashed lines: the low-redshift profile of lArnaud et al.l (120101). assuming se l f-simila r scaling 
to higher z. Dot-dashed lines: the mea n simul ated profile of lBattaglia et al.l (l2011bl ). Dotted 
lines: mean profile using Eqn.[8]for a (jSNBLl ). 
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Fig. 6. — Comp arison of the model to X-ray observations by ISun et al.l (120091 ) and 



Mantz et al.l (120101 ); for the latter, error bars have been omitted for clarity. The red line is 



the median and vertical bars enclose 68% of the model halos; at the high T end individual 
halos are shown as circles. 
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Fig, 



Planck Collaboration fc et al. 



Comparison of th e 
(boiltfl . 



Planck Collaboration fc et al. 



model to X-ray and 
Left: points with 



SZ observations by 
error bars are from 



( 12011b! ) . The red line is the median and vertical bars enclose 
68% of the model halos; at large Y individual halos are shown as circles. M g is measured in- 
side r 50 o and T in the radial range [0.15-0.75]r 50 o. All model halos in the light cone with z < 
0.5 are included. Right: SZ signal versus X-ray luminosity in t he band [0.1-2.4]keV, measured 



inside r 5oo- Squares with error bars are bi nned clusters from iPlanck Col 



fcOllaf ); crosses are individual clusters from lPlanck Collaboration fc et al. 
bars omitted for clarity 



aboration fc et al. 



( 12011b! ) with error 
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Fig. 8. — Projected cylindrical SZ decrement for halos hosting L RGs. Circles with error 
bars are from ACT data, with masses estimated from weak lensing (IHand et al.ll201ll ). The 
radio-quiet LRG sample is shown; masses are measured inside 200 times the mean density. 
The solid red line with vertical bars is the median and 68 % range of the f nt h model clusters. 
The dotted line is the median for the Sta ndard model of IBOV09I . and the dashed is for the 
nonthermal20 model of iTrac et al.l ( 120111 ). 



